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Abstract. We present a new staggered discretization of the Dirac operator. 
Doubling gives only a doublet of Dirac fermions which we propose to interpret 
as a physical (lepton or quark) doublet. If coupled with gauge fields, an (I+7 5 ) 
chiral interaction appears in a natural way. We define a generalization for 
curved background which does not require tetrad variables. The approach 
suggests a natural explanation for the three fermion families. 
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1. Introduction 

In this paper we present a new discretization of the Dirac equation. In com- 
parison with staggered fermions [8] it creates not four but only two flavours of 
Dirac fermions. This has been reached by placing not only different spin compo- 
nents, but also their real and imaginary parts into different nodes. These sixteen 
degrees of freedom (two fermions) can be understood, in some sense, as the result 
of "doubling" of a real scalar step of freedom ip(n) on the lattice. 

Moreover, these two fermions live on different sub-lattices (ip on "odd" nodes, 
ip e on "even" nodes), thus, we obtain also a single fermion (eight degrees of freedom) 
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on the lattice by omitting half of the lattice nodes. But, if the neutrino is a Dirac 
particle, all fermions appear in doublets of Dirac particles. In this context, the 
appearance of a fermion doublet in this discretization may be not a bug but a 
feature which allows to explain why Dirac particles appear in such doublets. 

In our approach the complex structure is an operator among others. Moreover, 
there is no natural complex structure, but, instead, a quaternionic structure. The 
choice of a complex structure corresponds to a choice of a direction. To avoid 
unnecessary symmetry breaking we propose to consider a vector field instead of a 
scalar field on the lattice. In this case, each component has its own natural preferred 
direction and therefore a natural complex structure. This gives in a natural way 
three generations of our doublet. 

Once a complex structure is defined, we can define the operator 7 5 on the grid. 
It exactly anticommutes with D. On the other hand, (7 s ) 2 ^ 1 but a shift by 2a 
in the preferred direction. The subspaces 7 s = ±1 have an uncertain definition but 
a clear physical interpretation: they define the subspaces ip w ±?/> e . An 1 + 7 13 
type chiral interaction with a gauge field appears in a straightforward way from 
the typical lattice gauge term A n<n + ai (<fi(ri) + 4>(n + Oj)): If 7 5 ~ 1 it becomes 
« Ai4>(n), but if 7 5 « — 1 it becomes almost zero. 

The discretization may be generalized to a curved metric background. The 
fermion doublet may be, in this context, described without any tetrad variables 
by the sixteen-dimensional exterior bundle of four-dimensional differential forms 
fi = X)fc=o^ fe - On this space exists a natural "Dirac operator" (in the sense of 
a square root of the harmonic operator) as well as a natural discretization for a 
general mesh. 

In appendix [X] we consider the connection of this concept with a generalization 
of Lorentz ether theory to gravity. In this context, the proposal to explain fermion 
families as natural degrees of freedom for ether crystal dislocations becomes really 
meaningful. Moreover, we define a natural SO(3) gauge field which measures mesh 
irregularity which we propose as a candidate for the SU (2) part of electroweak 
theory. 

There is a large number of proposals for chiral gauge theories on the lattice. They 
are mainly focussed around the Ginsparg- Wilson relation [5] or its generalizations 
[7], especially domain wall fermions [IB] . Neuberger's overlap operator [IT], propos- 
als by Fujikawa [4] and Chiu [3]. Because in our approach (^^esh) 2 1) there is 
no straightforward connection between our approach and these other approaches. 



2. A Real Representation of the Dirac Algebra 



We forget - for some time - about the complex structure. Instead of the usual 
representations with four complex fields, we use an eight-dimensional real represen- 
tation of the operators 7^ defined here by their linear combination with d^: 
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In the context of this representation, it seems also natural to define (by their 
linear combination with scalar parameters m,) the following operators /3 Z : 



(2) m^ = def 
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The following operator equation holds: 

(3) (7% - i l d, + m#) 2 = -□ + S^m^ 

This can be easily seen - this operator iterates three times, in each coordinate 
direction, the same trickQ 



(4) 
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As follows immediately, the 7 M define a representation of the Dirac matrices, and 
the matrices (5 l fulfil the following anticommutation relations: 

(5) I3 l (i j + (3 j f3 l = 5 l] 

and anticommute with all 7 M : 



(6) /3V + 7^/3 4 = 

It is also easy to see (and to generalize to arbitrary dimension) that 

(7) 7 °(7 1 /? 1 )(7 2 /3 2 )(7 3 /? 3 ) = l. 

The "classical" operator 7 s = — 17°7 1 7 2 7 3 depends on the complex structure, 
which is not yet defined. The natural replacement which does not depend on it 
the expression 7°7 1 7 2 7 3 - we denote with l: 



This observation also suggests how to iterate this construction to arbitrary dimension. 
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(8) i = def 7V7V = /3 1 /? 2 /? 3 ^+7^ = W 2 = -l 

3. Discretization of the Dirac Equation 

This representation is appropriate for a discretization of the Dirac equation on 
a regular hyper-cubic lattice. It can be obtained in a quite simple way: We start 
with a naive central difference approximation 

(9) dii/)(ri) — * — Upin + cti) — ip(n — en)). 

This naive discretization leads to the problem of "fermion doubling" . The con- 
tinuous limit of this set of discrete equations gives not only the original Dirac 
equation, but also additional, highly oscillating components, the "doublers". In 
classical computations such doublers may be often ignored, but in quantum com- 
putations, where the number of degrees of freedom is important (Pauli principle) 
this is no longer possible. We obtain in each direction a factor to, thus, 2 4 = 16 
doublers. Fortunately, eight pairs of doublers decouple in a really simple way: It is 
sufficient to hold only one of the eight real components i\) a per node. On the three- 
dimensional reference cube (e\, E2E3)) £% € {0, 1} we obtain the following locations 
for the eight components: 

t[)° located at (0, 0, 0), V 4 located at (0, 0, 1); 

, . V 1 located at (1,0,0), V 5 located at (1,0,1); 

^ ' V 2 located at (0, 1,0), V 6 located at (0, 1, 1); 

V> 3 located at (1, 1,0), V 7 located at (1, 1, 1); 

What remains are sixteen degrees of freedom (eight degrees of freedom on two 
time steps which we need because of our use of central differences in time) which 
corresponds to a doublet of Dirac fermionfl Note that our discretization may be 
interpreted as a way to discretize the d'Alembert equation for a single scalar step 
of freedom ip(n) with central differences, which gives 2 4 = 16 doublers. 

Now, the last doublet decouples too, but in a slightly less trivial way: We can 
distinguish "even" and "odd" nodes on the full space-time lattice. The central 
difference equations on even (odd) nodes connects only values on odd (even) nodes. 
Thus, we obtain two fermions ip e and ip on even resp. odd nodes. On the four- 
dimensional reference cube (so, £1, £2£a), £j € {0, 1} we have 

ip e ° located at (0, 0, 0, 0), Ve 4 located at (1, 0, 0, 1) 

, , tpe 1 located at (1, 1, 0, 0), ip e 5 located at (0, 1, 0, 1) 

^ > Ve 2 located at (1, 0, 1, 0), ?p e 6 located at (0, 0, 1, 1) 

ip e 3 located at (0, 1, 1, 0), ip e 7 located at (1, 1, 1, 1' 



2 This is a variant of a well-known approach to solve the "fermion doubling" problem - so- 
called staggered fermions [5]. This approach reduces the doublers only by a factor four. Because 
we ignore the complex structure of the standard representation, we are free to place "real" and 
"imaginary" part of the complex fields into different nodes. This gives the additional reduction 
by factor two. 
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located at (1, 0, 0, 0), Vo 4 located at (0, 0, 0, 1); 

, . Vo 1 located at (0, 1, 0, 0), tp 5 located at (1, 1, 0, 1); 

^ ' Vo 2 located at (0, 0, 1, 0), Vo 6 located at (1, 0, 1, 1); 

tpo 3 located at (1, 1, 1, 0), ip 7 located at (0, 1, 1, 1); 

But, instead of removing one sub-mesh to describe a single Dirac fermion, we 
propose to accept above doublers as a way to describe a physically meaningful 
flavour doublet. Remarkably, if the neutrino is a standard Dirac particle, then all 
fcrmions of the standard model appear in doublets. The appearance of a fermion 
doublet in our approach may be, therefore, not a bug but a feature which allows to 
explain the existence of these doublets. 

4. Complex Structures 

To connect the Dirac particle with a gauge field we need a complex structure - 
a multiplication with i. In our real representation this is an operator. It's basic 
property is i 2 = — 1. To be compatible in the usual way with the Dirac equation, 
we also need [7^,1] = 0. This does not define the complex structure uniquely. We 
have several interesting candidates for a complex structure: 



(13) i ^P 1 ^ 2 =l(3 3 

(14) j = (3 2 f3 3 =i$ x 

(15) k = {3 3 [3 1 =l(3 2 
which together define a quaternionic structured 



(16) ij — —ji — — fc; jk = — kj = —i; ki = —ik = —j; i 2 — j 2 = k 2 = —1 

For each candidate for a complex structure, we obtain an own operator 7 s —def 
—it. Especially for i — l(3 3 we obtain 7 s = — l(3 3 l = f3 3 . 

4.1. Fermion Families and Lattice Distortions. Thus, it seems that to fix the 
complex structure we somehow have to break spatial symmetry. But there is a 
simple way out of this. Instead of one scalar step of freedom ip(n) on the lattice, we 
can consider a vector field - thus, three components f l {n). Now, each component 
has a natural "preferred direction" and, therefore, a natural complex structure. 

Moreover, a vector on a lattice is a quite natural step of freedom. It is, for 
example, the natural way to describe lattice distortions with a shift vector field 
u l (n). 

On the other hand, in the standard model we have three families - three copies 
of each fermion. This suggests to explain on the kinematic level the three fermion 
families using the hypothesis that the fundamental degrees of freedom of the "the- 
ory of everything" are three-dimensional vector fields u l (n). This idea becomes 
really meaningful only in an ether-theoretical context, where relativistic symmetry 
is considered as not fundamental but an effective large scale symmetry, and dis- 
tortions of an "ether crystal" are natural degrees of freedom. We consider such an 
approach in more detail in appendix [A] 



The classical representation ij = k can be obtained using reverse signs for i,j, k, but we prefer 
this sign convention because it gives 7 s = /3 3 . 
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For our symmetry problem - the choice of the complex structure - other solutions 
are possible. We can use them all in expressions like 



(17) iA 1 + jA 2 + kA 3 

which allows to connect with an SU{2) or SO (3) gauge field. We can also try 
to connect colors with spatial directions to define the preferred complex structure 
differently for each colored quark. Last not least, it may be really broken (possibly 
the U{\) part of electroweak interaction). These different solutions possibly even 
do not contradict each other. 

But it seems important to notice that the choice of the complex structure is not 
natural, requires a symmetry breaking of the quaternionic structure. This point 
may be important for the understanding of symmetry breaking in the standard 
model. 



5. Chiral Symmetry on the Lattice 

Let's assume now that one of the complex structures, namely i — t/3 3 , has been 
chosen. To understand chiral symmetry we have to define the operator j 5 = (3 3 
on the lattice. Note that it cannot be anymore a pointwise operator as for Wilson 
fermions and staggered fermions - it connects components which are located in 
different points. Thus, the original, continuous definition is not enough to define 
7 5 on the lattice. But there is one very natural choice: the shift operator in z- 
direction. Modulo 2, on the reference cube, it is indeed j 5 . And, most important, 
it anticommutes with the Dirac operator even on the lattice: 



(18) 7 5 L> + L> 7 5 =0 

But, different from the j 5 in continuous theory, it's square is no longer 1, but a 
shift in z-direction. This is a very nice and beautiful way to break the exact chiral 
symmetry (7 5 ) 2 = 1. Note that 7 5 exchanges the even and odd fermion: 



(19) 7 5 -0 e (n) = ipo(n - a z ) -f 5 ip (n) = ip e (n - a z ) 

This gives the subspaces ip± = ^(1 ± 7 5 )-0 a quite obvious physical meaning: 
they are the subspaces defined by ip e w ±ip a . On the other hand, because tp e and 
ipo are located on different sub-meshs, we do not have a natural symmetric and 
exact definition of the projectors. 

Now, the most beautiful surprise we find thinking about the connection to gauge 
fields. As usual, we describe gauge degrees of freedom as located on edges. Now, 
a quite natural, "naive" symmetric lattice interpolation for interaction terms Atp 
between fermions and the gauge field is 



2 

Thus, we obtain a dependency on terms of type ip(n + at) + ip(n) in a natural 
way. Now, these terms are w ip(n) for ip + but w for ip_ - thus, our interaction 
terms with gauge fields have exactly the (1 + 7 5 ) form of interaction we need in 
chiral gauge theory! Thus, in this discretization we have not only a possibility to 
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describe chiral gauge fields, but they appear almost by themself, in a quite natural 
way. 

This effect disappears if we have several components and the gauge field connects 
them. Thus, our (1 +7 5 ) type interaction appears only inside the doublet. Outside 
the doublets we obtain a vector-like interaction - as it should be in chromodynamics. 

6. The Dirac Operator on the de Rham Complex 

The discretization may be generalized to the case of a general metric background 
in a quite simple way - as the exterior bundle or the de Rham complex. The de 
Rham complex and it's Dirac operator is well-known, but has the "wrong" dimen- 
sion for a single fermion: in four dimension, we have 2 4 = 16 components. Instead, 
fcrmions are usually described on a curved background using tetrad variables for 
the metric (see, for example, pQ). 

In our approach we have to describe a doublet of Dirac fermions, and in this 
context these sixteen components of the exterior bundle is exactly what we need. 
In this sense, the "fermion doubling" problem has a natural analogon in continuous 
theory on a curved background: only a pair of fermions may be described on a 
curved metric background g^ v without tetrad variables. 

6.1. Hodge Theory. Let's remember the basic formulas for the Dirac operator 
in the exterior bundle (see, for example, [14]). The exterior bundle or de Rham 
complex il — J2k=o ^ k cons ists skew-symmetric tensor fields of type (0, k), < k < 
n which are usually written as differential forms 

(21) ij) = ipi 1 ...i k dx il A • • • A dx ik G Vt k 

The exterior bundle Q has dimension 2™ in the n-dimensional space. The most 
important operation on fl is the external derivative d : Sl fc — > £l fc+1 defined by 

( 22 ) (^^-Eif- 1 )^...:,^, 

9=1 

where i q denotes that the index i q has been omitted. It's main property is d 2 = 0. 
In the presence of a metric, we have also the important ^-operator £l k — ► f2"~ fc : 

(23) (*V0w..<„ = ^n.., n 9 lin ■■■g lkJk Tpj 1 ...3 k 

with * 2 = (— l) fc (™~ fc )sgn(g). This allows to define a global inner product by 



(24) (<^) = / MH0 = / V- A (*</>) 
It turns out that the adjoint operator of d* of d is 

(25) d* = (-l) rn+n+1 * d* 

In this general context we can define the Laplace operator as 



(26) 



A = dd* + d*d 



s 



ILJA SCHMELZER 



Then, the Dirac operator (as it's square root) can be defined as 



(27) 



D = d + d*. 



Indeed, we have d 2 = as well as (cT) 2 = 0. Sometimes the Z2 graduation is 
also useful: eif> = (—l) k ip if ip G Q k ■ The subspaces e = 1 and e = — 1 have equal 
dimension 2"~ 1 . 

6.2. Discretization on the de Rham Complex. If we want to consider the 
approximation of some continuous object on a general mesh, the exterior bundle 
is a very natural object. Indeed, k-forms may be integrated over k-dimensional 
surfaces. Thus, if a general mesh (a cell complex) is given, a k-form defines a 
function on the k-dimensional cells of the mesh in a very natural way. For each cell 
we can define it's "characteristic form" on the mesh as the equivalence class of forms 
with x l { c j) = $j an d decompose the general form on the mesh as /(c) = fiX l {c). 
The external derivative defines in a similar natural way a derivative for functions 
on the mesh, with the same most important exact property d 2 = 0. 

For the definition of the *-operator we need a metric and a dual mesh. The 
metric defines in a natural way for every cell Cj it's area a; = a(ci) > 0. In the 
Euclidean case and a triangulation, these values depend on each other, but in the 
general case they may be considered as independent variables which define the 
metric. In the following we consider them as given and defining the metric. 

The dual mesh is a quite natural mesh which can be defined for every mesh: for 
each k-dimensional cell Ci of the original mesh it has a (n — fc)-dimensional dual cell 
cj which intersects the original cell in a single point orthogonally and with positive 
intersection index. Now, the metric defines as well a scalar value a, for each cell of 
the dual mesh £i - the area of the dual cell. Now, the *-operator on the mesh may 
be defined as 



Note that the dual of the dual mesh has the same cells, but in some dimensions 
with different orientations. Therefore, for * 2 there appears the factor (— i)' c (™- fe ) 
as in the continuous case. 

Note also that the discrete exterior bundle is connected in a simple way with 
a natural "regular refinement" step of a mesh: To refine a mesh in a regular way, 
we can put one node in the center of each cell, and connect nodes of neighbour 
dimension cells with edges if one of the cell is on the boundary of the other. This 
leads to a mesh where each cell has the topological structure of a unit hypercube. 

For a simple hypercubic lattice on a flat background, this discretization scheme 
reduces to our proposal above. 



The approach described here may help to solve some outstanding problems of 
chiral gauge theory. First, there is the general problem how to formulate a chiral 
gauge theory on the lattice, known as the "regularization problem" : any consis- 
tent regularization that preserves the gauge symmetry must refer to the fermion 
representation. 12J. In our case, we describe chiral gauge theory only in a very 



(28) 



(*/)(c) = -rfix'ic) 



7. Conclusion 
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special, two-dimensional representation, therefore it is not in conflict with this re- 
quirement. And it gives a representation we need in the standard model to describe 
the SU (2) part of the electroweak interaction. Moreover, if neutrinos are Dirac par- 
ticles, we don't need any other representations - only their combinations (one lepton 
and three quark doublets). Of course, a lot of things have to be done to obtain a 
complete quantum chiral gauge theory in this way. The problem how to handle the 
U(l) part of electroweak interaction, which needs the incorporation of all fermion 
doublets, is completely open. 

Second, it may be helpful for a solution of the fermion doubling problem (see 
[6] for an introduction, [7,9,12,17] [9] [12j [17] for various approaches to solve 
this problem in chiral gauge theory). It reduces the number of "doubled flavours" 
from four in the staggered fermion approach (which seem unphysical) to two which 
already allow a meaningful physical interpretation as lepton or quark doublets. 
The nontrivial character of the complex structure may prevent the application of 
standard no-go-theorems like the famous Nielson-Ninomiya [10] theorem. 

We also obtain an interesting representation of fermions on a general metric 
background which does not need tetrad or triad variables for the gravitational field. 
Instead, the gravitational field may be described as usual by the metric tensor g^. 

8. Acknowledgements 
Thanks to M. Nobes, L. Bourhis and S. Carlip for interesting discussions. 

Appendix A. Connection to the General Lorentz Ether 

The approach to fermions presented here has been developed to solve a particular 
problem connected with the description of fermion fields General Lorentz Ether 
Theory (GLET), which we want to describe here shortly. 

This theory has been proposed by the author in [15] as a generalization of the 
classical Lorentz ether to gravity. It is a metric theory of gravity, where the "effec- 
tive" metric g^ v is defined by a variant of the ADM decomposition from classical 
condensed matter (ether) variables: density p = g°°\/— g > 0, velocity v % = g 0l /g 00 , 
and a three-dimensional pressure tensor p lJ y/p = —g lJ ^fg. The continuity and Euler 
equations of classical continuum mechanics 

(29) dtp + diipv*) = 

(30) d t (pv j ) + d l (pv i v j +p ij ) = 0. 

transforms into the harmonic condition d^g^yf—g = for the effective metric. 
The Lagrangian of GLET consists of the classical GR Lagrangian and some non- 
covariant term which describes the dependence on the Newtonian background: 

(31) L = L GR + L matt er{g^^ m ) - (87rG)- 1 (T«? 00 - 3S ijg ^)V=g, 

Therefore, the theory has a well-defined GR limit T, 3 — > which makes the the- 
ory viable. This Lagrangian has been derived from first principles which seem quite 
natural for condensed matter theories in [15] . Especially the covariant character of 
its matter part, that means the Einstein equivalence principle, follows. 

It has been noted [5] that this derivation does not extend (at least not in any 
obvious way) to tetrad or triad variables. Because the standard way to describe 
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fermions on a curved background uses tetrad (or, in the context of an ADM decom- 
position, triad) variables, this makes the incorporation of fermions into this theory 
problematic. One way to solve this problem is to find a way to describe fermions 
without tetrad or triad variables, what has been done in our approach: Here we 
define a pair of Dirac fermions as classical tensor fields. 

It should also be noted that the regularization problem of chiral gauge theory 
as well as the fermion doubling problem become important in ether theory if we 
want to use an "atomic ether" concept to regularize the theory. An "atomic ether" 
regularization, which interprets the density p as a particle density and velocity v % 
as their average velocity is a quite natural way to quantize an ether theory. But 
in this case we need some lattice discretization for all parts of the standard model 

In the context of the ether approach, an extremely simple toy model would be a 
cubic ether crystal. Distortions of this crystal would be defined by three real values 
per node. This toy model exactly fits into the lattice discretization described here, 
which gives three generations of a pair of fermions. Thus, already an extremely sim- 
ple toy model seems to be able to explain the three fermion generations. Moreover 
it suggests that fermions appear in pairs, with a natural possibility for 1 + 7 s type 
gauge interactions inside these pairs. These are already nontrivial and important 
features of the standard model. Especially if neutrinos appear to be Dirac particles, 
this scheme holds for all fermions of the standard model. 

Appendix B. Compatibility with ADM Decomposition 

For the compatibility of the approach described here with GLET it is the com- 
patibility with the ADM decomposition is crucial and worth to be considered. Using 
"comoving" spatial coordinates (which remain constant along the velocity field v l ), 
the harmonic operator of the metric g^ v reduces to 

(32) [ty = -(p#-A)V 

where A is the standard three-dimensional (harmonic) Laplace operator of the 
metric pij . Thus, we ADM decomposition allows to reduce the harmonic equation 
to a "three-and-one-half-dimcnsional" equation 

(33) y/pdttj) = \/A>. 

Now it seems natural to use a mesh which is compatible with the ADM de- 
composition. We start with an arbitrary mesh on the spatial slice at a given time 
to- On this three-dimensional mesh we can define the Dirac operator following the 
scheme described above. Then we can extend this mesh into a four-dimensional 
mesh. First, we extend the mesh nodes at to following the comoving spatial coor- 
dinates into comoving world-lines. For the preferred time we use a global, regular 
lattice and simple central differences. Such a choice of a lattice and discretization 
is in natural agreement with a crystallic ether model. 

Appendix C. A Natural SO (3) Gauge Field for a General 

Three-Dimensional Mesh 

Another reason why we think that the focus of three spatial dimensions will be 
helpful is that a general three-dimensional mesh defines in a simple and natural 
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way an SO(3) gauge field. Because at the Lie algebra level so(3) = su(2) this may 
be a candidate for the su(2) gauge field of the electroweak interaction. 

To understand this gauge connection, remember that to define a gauge connec- 
tion it is sufficient to define its integral over closed loops. Thus, let's define the 
integral over a closed loop starting from a general point (that means, inside a cell) 
and of general type (thus, the loop intersects only planar cells of the mesh). 

Now, assume also that we have made a simple regular refinement step, therefore, 
the mesh is a cubic mesh. For each cubic cell, we fix some standard map into the 
regular cube. Then we start with a standard reper (as it comes from the standard 
cube) and use trivial parallel transport of the reper (as defined on the standard 
cube) inside the cells up to the boundary. On the boundary, one vector of the reper 
is always orthogonal to the boundary. If we come back to the original node, the 
transferred reper is also, by construction, a reper in this cube, its directions are 
parallel to the original reper. Nonetheless, it is in general a rotated reper. This 
rotation defines the gauge connection. It measures the irregularity of the mesh: For 
a regular cubic mesh it is trivial. 

Now, the existence of a gauge field which is in such a natural way connected with 
mesh irregularity does not look like an accident of nature. Our hypothesis is that 
it is (at least connected in some way with) the SU(2) gauge field of the electroweak 
interaction. The details of such a connection are subject of future research. 
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